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ABSTRACT 

We study the effect of corotation resonance on the inertial-acoustic oscillations (p- 
modes) of black-hole accretion discs. Previous works have shown that for barotropic 
flows (where the pressure depends only on the density), wave absorption at the coro- 
tation resonance can lead to mode growth when the disc vortensity, C = k 2 /(2J7E) 
(where Q, n, S are the rotation rate, radial epicyclic frequency and surface density 
of the disc, respectively), has a positive gradient at the corotation radius. Here we 
generalize the analysis of the corotation resonance effect to non-barotropic fluids. We 
show that the mode instability criterion is modified by the finite radial Brunt- Vasala 
frequency of the disc. We derive an analytic expression for the reflectivity when a 
density wave impinges upon the corotation barrier, and calculate the frequencies and 
growth rates of global p-modes for disc models with various cv-viscosity parameteri- 
zations. We find that for disc fluids with constant adiabatic index T, super-reflection 
and mode growth depend on the gradient of the effective vortensity, <^ c g = 
(where S = P/S r measures the entropy): when d£ e g/dr > at the corotation radius, 
wave absorption leads to amplification of the p-mode. Our calculations show that the 
lowest-order p-modes with azimuthal wave number m = 2, 3, 4, . . . have the largest 
growth rates, with the frequencies approximately in (but distinct from) the 2:3:4... 
commensurate ratios. We discuss the implications of our results for the high-frequency 
quasi-periodic oscillations observed in accreting black-hole systems. 

Key words: accretion, accretion discs - hydrodynamics - waves — black hole physics 
- X-rays: binaries 



1 INTRODUCTION 

High frequency quasi-periodic oscillations (QPOs) in X-ray binary systems have been observed for a number of years and 
may provide an important tool for studying the strong gravitational fields of black holes (see Remillard & McClintock 2006). 
However, the physical mechanisms that generate such X-ray variability remain unclear. One of the most appealing models 
for the source of QPOs is the relativistic diskoseismic oscillation model, where general relativistic effects produce trapped 
oscillation modes at the inner region of an accretion disc (e.g., Kato & Fukue 1980; Okazaki et al. 1987; Nowak & Wagoner 
1991; see Wagoner 1999 and Kato 2001 for reviews). Other related works on black-hole diskoseismology, such as possible mode 
excitation and damping (e.g., Ortega-Rodriguez & Wagoner 2000; Li, Goodman & Narayan 2003; Kato 2003,2008; Tagger & 
Varniere 2006; Ferreira & Ogilvie 2009; Tsang & Lai 2009a), the effects of disc magnetic fields (e.g., Tagger & Pellat 1999; Fu 
& Lai 2009) and numerical simulations (e.g., Arras et al. 2006; Reynolds & Miller 2008; O'Neill, Reynolds & Miller 2009), as 
well as other ideas for high-frequecy QPOs, such as non-linear resonances (e.g. Abramowicz & Kluzniak 1999; Horak & Karas 
2006; Rebusco 2008) and boundary layer oscillations (e.g. Li & Narayan 2004; Tsang & Lai 2009b), are reviewed in section 1 
of Lai & Tsang (2009). 

* Email: dtsang@astro.cornell.edu; dong@astro.cornell.edu 
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In Lai & Tsang (2009), we studied the global corotational instability of non-axisymmetric p- modes (also called intertial- 
acoustic modes) trapped in the inner-most region of the accretion disc around a black hole. These modes do not have nodes 
in the vertical direction, and were shown to be amplified by the effect of wave absorption at corotation resonance. Near the 
black hole the radial epicyclic frequency k reaches a maximum and goes to zero at the innermost stable circular orbit (ISCO). 
This causes a non-monotonic behavior in the fluid vortensity, £ = K 2 /(2Ef2), such that d(/dr > inside the radius where £ 
peaks. It can be shown that the sign of the corotational wave absorption depends on the sign of the vortensity gradient dQ/dr 
(Tsang & Lai 2008; see Goldreich & Tremaine 1979). Thus p-modes with positive vortensity gradient at the corotation radius 
can be overstable due to corotational wave absorption. Tagger & Pellat (2002) and Tagger & Varniere (2006) showed that the 
global p-mode intability can be enhanced when the disc is threaded by a strong (of order equipartion) , large-scale poloidal 
magnetic field. 

Our previous study (Lai & Tsang 2009) and much of the related work on disc dynamics have assumed barotropic 
flows for the disc (i.e. the pressure depends only on density). This assumption provides convenient simplification, but may 
miss important effects of the disc dynamics. For example, Lovelace et al. (1999) and Li et al. (2000) studied the adiabatic 
perturbations for waves trapped by a disk entropy radial profile that has a localized maximum, leading to the so-called 
Rossby-wave instability. They showed that in such a case the key parameter determining the effect of the corotation is no 
longer the gradient of the vortensity, but rather the slope of a modified effective vortensity, £ e ff = C,/S 2 ^ T , where S = P/E r is 
defined as the entropy and V is the 2-dimensional adiabatic index (assumed to be constant). As another example, Baruteau 
& Masset (2008) showed that the corotation torque of a protoplanetary disc on a planet can be significantly different for 
barotropic and non-barotropic fluids. 

In this paper we study the global corotational instability of p-modes in accretion discs around black holes, generalizing 
our pevious works (Tsang & Lai 2008; Lai & Tsang 2009) to include non-barotropic effects. In section 2 we develop the basic 
equations of adiabatic perturbations for generic accretion discs. In section 3 we analyze the effect of the corotation resonance, 
including a careful treatment of both the first and second-order singularities of the resonance; we derive a WKB expression 
for the reflectivity due to the corotation barrier and show that super-reflection can be achieved under certain conditions. In 
section 4 we consider black hole disc models parametrized by a-viscosity, and calculate the global disc p-mode frequencies 
and growth rates. Finally, in section 5 we discuss the implications of this work for models of high-frequency QPOs. 



2 BASIC EQUATIONS 

We begin by considering the basic fluid equations of a 2-dimensional disc. The continuity and momentum equations read: 
aS + u-VE + SV-u = 0, (fa) 
ftu+(u-V)u = -IvP-V<E>, (lb) 

where P(r) = j pdz is the vertically integrated pressure, and E(r) = f pdz is the surface density and we adopt the Pacyznski- 
Wiita psuedo-Newtonian potential $ = GM/(r — 2r g ) with r g = GM/c 2 . Assuming that the background flow has u = rQcj>, 
and that the Eulerian perturbations <SE, SP, and Su = 5u r f + 8u^(f>, have the form exp(im^) — iwt), we find the linear 
perturbation equations: 

f d Z7TL 

-iu5'E + - — (Y;r8ur) + — E<5u<* = 0, (2a) 
r or r 

-iCoSur - 2Q8 U4 , = --—8P+——, (2b) 
E Or E^ Or 

._ c k 2 . im 5P .„ . 

-lUlSu^ + — 5u r = — -|T , (2c) 

where c 2 — 9P/9E is the adiabatic sound speed, w = u — mQ, and k is the radial epicyclic (angular) frequency. 

Tsang & Lai (2008) and Lai & Tsang (2009) assumed the disc fluid is barotropic such that P — P(S). Here we consider 
adiabatic perturbations of a general non-barotropic disc. The Lagrangian density perturbation AE and pressure perturbation 
AP are related by 

AS = lAP . (3) 

This gives 

where £ r = iSu r /uj is the Lagrangian displacement in the r direction, and N r is the radial Brunt- Vaisala frequency as given 
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by 



N 2 = 1 



dP 
£ 2 V ~dr 



dE 
IP 



(5) 



Combining the above with the linearized perturbation equations (|2a[) - (|2c|) . and eliminating Su^, we obtain two coupled first 
order ODEs appropriate for numerical integration 

Sh' 



iSu' 



P> 

„,2 



cur J 



Sh 



Sh H — ^-idu r 



EN 2 



+ 



(lnr£)' 



iSu r , 



where Sh = <5P/E is the enthalpy perturbation, " ' " denotes d/dr and 



D s = K 2 — LU 2 + N 2 



(6a) 
(6b) 

(7) 



Eliminating 8u r we arrive at the second order differential equation for Sh, 







where 



<f_ 5h _± 
dr 2 dr 



hi 



rE J dr 



+ 



dr 



Sh- 
1 

17 



2 ' r .2 

j__d_ 

L s dr 



2mO d 



In 



r<I> dr 
D 



hi 



ns 
#7 



Sh 



rE 

1 
17 



4mft m 2 iV 2 
+ t 2 - 

urLs r 2 w 2 



EN 2 



Sh 



dP/dr ■ 

It is convenient to eliminate the term proportional to Sh' in eq. © by defining 



A 2 ^ 
rE 



and 



S/i 



(8) 



(9) 



(10) 



which allows us to rewrite 



<9r 2 



m D s 2mf2 d 



as a wave equation 



J_ d_ 

If + dr 



ro; dr 



In- 



1 
IT 



J__d 

Ls dr 



hi 



rE 



ill 

dr 2 A j 

4mfi m 2 AL 2 



ujrLs 



(11) 



Equation (|llf) forms the basis of our analysis in section 3. If the adiabatic index T = <91nP/<91n£ = c 2 E/P is constant 
one can define the "entropy" , 



Then 



and 



LS 1 



S = P/E 1 

1 dP dlnS 

" _ fE"d7 dr 

1 d InS 



T dr 



(12) 
(13) 
(14) 



and eq. ((8)l reduces to eq. (10) in Lovelace et al (1999). When N r — > (thus L a 1 — > 0), the terms on the second line of eq. ([8} 
vanish and we recover the second order perturbation equation for barotropic flows (Goldreich & Tremaine 1979; Tsang & Lai 
2008). 



3 REFLECTION OF THE COROTATION BARRIER 

Away from the corotation resonance (where to = 0) region, eq. (|lip yields local WKB wave solution Sh oc exp(i J k r dr), with 
D s /c 2 ~ —A: 2 , or 

uj 2 ~ k 2 + N 2 + k 2 c 2 . (15) 

Since typically N 2 <^ c 2 /r 2 <C k 2 ~ £2 2 (for thin discs), this is the standard dispersion relation for spiral density waves. 
The inner/outer Lindblad resonances (I/OLR) are defined by u = ±-\/k 2 + N 2 ~ ±k. Waves can propagate inside the ILR 
(r < til) or outside the OLR (r > roL). between t*il an( l r OL lies the corotation barrier. 

In this section, we derive the expression for the (complex) reflection coefficient for waves incident upon the corotation 
barrier and deduce the condition for super-reflection. Our analysis generalizes that given in Tsang & Lai (2008), which assumed 
barotropic fluids. 
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3.1 Analytical Calculation of the Reflectivity 

Near the corotation resonance r = r c where Co = 0, we can rewrite eq. as 



dr 2 



q\dr fiE L s J r ~ R c q 2 Q 2 (r ~ R c 



where q = — (dlnS7/dlnr) c , R c = r c — i^^, and k e g given by 



2 D, 



d 2 1 



d 



-r + ^-A—- + -^ + - — -__ Mn-i 



1 d 



Ls dr \ rE 



»7 = , 



(16) 



(17) 



r^ cj dr 2 A L s dr 
is the effective radial wave number without the terms singular at the corotation. We have introduced a small imaginary part 
to the wave frequency so that uj = ui r + iut (with uji > 0). Defining 

2k eS r c uji 



2fc cff dr, 



and e 



we have 



dx 



1 



+ 



1 2 

4 "A* 



V> = 



4 ' x + it ' (x + ie) 2 

which we recognize as the Whittaker differential equation (Abramowitz & Stegun 1964). In eq. (|19fl we have defined 



d_ 

qn \ dr 



InC- 



2 



ill 



AN 2 
q 2 Q. 2 



2 \ 1/2 



where 



2ftE 



(18) 
(19) 

(20) 
(21) 

(22) 



is the vortensity for the background flow. When the adiabatic index F = constant, we can use eq. (|14[1 for Lg 1 , and define 
the effective vortensity f e ff so that 



c s d 
qn dr 



with 



Ceff ; 



(23) 



(24) 



2fiES' 2 / r ' 

Typically « ~ £1 (away from the ISCO) and |JV r | ^ c s /r <C (for thin discs), and we have in order of magnitude \u\ ~ c s /(rQ) 
and ~7V 2 /c 2 ^c 2 /(rf7) 2 . 

Equation (|19[1 is solved by the Whittaker functions with indices v and fi. The two linearly independent functions of 
z = x + ie convenient for construction of connection coefficients are 



and 



'W- v ^(ze- tlT ) + -T W„, M (z), 



(25) 



where To is the stokes multiplier (defined below), and z is defined such that arg(z) ranges from to it. The resulting connection 
formulae (Tsang & Lai 2008) are: 



5h- 



5h A 



I fccff 

_A 

A 



exp 



e^exp (+/; c fc cff dr) + . 

exp ( + k e g dr^ 
^e'^exp (+/; c fc off dr) 



It. 



™ exp (-/;<= fc cff dr) 



for r»r c 
for r <C 



where the Stokes multipliers (Heading 1962) are given by 

To 



for r>r c 

' exp (- f r " fc cff dr) for r < r c 
2 7 rie l2 '"' 



(26) 
(27) 

(28) 



The connection formulae (|26|) - (|27[ ) here are the same as eqs. (42)- (43) of Tsang & Lai (2008), the differences lie in 
the expressions for v [eq. (|20[) ]. To and T\ [eq. (|28[) ]. For barotropic fluids, Ls — > oo and — > (and thus fi — > 1/2), our 
expressions reduce to those given in Tsang & Lai (2008). Using the connection formulae ()26p - (|27p and the connection formulae 
for the Lindblad resonances given by eqs. (32)-(35) in Tsang & Lai (2008), we can obtain the reflection and transmission 
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coefficients for waves incident upon the corotation barrier for r < til: 



1Z : 



l+i e -;2™ e -2e„ ^ + i T()Tl ) + ^ Tie -i2^ e -2e Ua _ | Toe -2e IIb 
1 - I e- i2 ™e- 2@ n (1 + iT Ti) - frie- i2 ™e- 2e na - f T e- 2e nb 



r = 



1 - 1 e-^^e-aen (i + |t Ti) - iTie-^e- 



where 



9t 



/ r ° L /^dr, 6„ a = r 

r T T. J TTT. 



iToe- 2e nb 



— k 2 „dr, and 



e nb = 



For |^| <C 1 and /i — || <C 1 eq. (|29l) can be simplified to 



ft = 



e e " + Ie" ei 



-20 



i e -e : 



+ 



IIb + i - %e 



1 „-2©r 



-2H 



( e ©n - i e - e n 



-7TI/ + 



lib _ 2@llb 

4 e 



/"■ol / 

/ \J-k 2 eS dr. 

7T( M -i)+0[^,( M -i) 2 ] 



For 0nb S> 1 and 8n a 3> 1 this further reduces to 



R - I - I v + fj, - I ) 7re- 2eiIa 



Thus, super-reflection occurs when ^ + /i + i~^>0 [since |/Lt — |J is much smaller than see eqs ((20)) - (f2T|) ] 



(29) 
(30) 

(31) 

(32) 
(33) 



3.2 Numerical Calculation of Reflectivity 

We can also calculate the reflectivity numerically by integrating eqs. (|6a|) - (|6bjl . To this end, we assume an outgoing wave at 
some radius r out 3> i~ol, motivated by the wave equation 



<5/i oc 

where fc r is the full radial wave-number given by 
a D a 2mQ, d (, fiE 



A 



exp 



k r dr + i- 



^ k 2 r = — + ^ 
r r 2 c? 



d 2 1 



This gives the outer boundary condition at r = r D 



dr 
1 d 



J__d_ 

is dr 



rE J d)rl/s 



4mQ, m 2 N 2 
H 11 

1 — o 



5/i'(r out ) = ( ik r + — In A — -— lnfc, 
1 dr 2 dr 



<5/i(r 



At some inner radius n n <C riL, the solution takes the form 



Sh 



A 



exp 



k r dr + i- 



+ 72. exp 



fc r dr 



and the reflection coefficient can be obtained from 



\n\ = 



■In A 



l_d_ 

2 dr 



In k r — ifcr) <5/i — Sh' 



-feln A -±A mfcr + ikr ) 5h-6h> 



(34) 



(35) 



(36) 



(37) 



(38) 



Fig-Dgives some examples of the reflectivity for a simple power-law disc model. We parameterize the relevant disc profiles 



by assuming E oc r p , N r = N c (r/r c ) and Ls = Ls c (r/r c 



and we fix the sound speed to c a /(rQ) =0.1. Note that 



for non-barotropic fluids P 7^ Ec 2 , thus we allow N 2 and Lg 1 = T,N 2 / P' to vary independently. In the examples depicted in 
Fig. [TJ we fix N r (r) and Ls(r), but vary the density index, p, to change the critical parameter v [see eq. (|20|) ]. In general, as 
seen from eqs. (|6al) - (|6b[) or eq. (|11[1 . the result depends on 7V C , Ls c , b, and 7, but we find that the dependence on b and 7 to 
be rather weak. In agreement with the analytic expression in Section 3.1 [eq. I|33[)]. we find that the super-reflection > 1) 
is achieved when 1/ + n — |~^>0. 



4 CALCULATION OF GLOBAL OVERSTABLE P-MODES IN BLACK HOLE ACCRETION DISCS 

The result of Section 3 shows that when v + fi — \ > 0, waves impinging upon the corotation barrier are super-reflected. 
Supposing there exists a reflecting boundary at the inner disc radius n n = nsco, normal p-modes can be produced, with 
waves trapped between ri n and hl. In the WKB approximation the mode growth rate uji is directly related to the reflectivity 
\TZ\ (Tsang & Lai 2008) 

~<M)\Cd5k>*r (39) 

where ui r — Re(cj) — mfi. More accurate calculation of the p-mode frequency u> = ui r + iuji requires solving the complex 
eigenvalue problem based on eqs. (|6a|) - ()6b|l . 
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Figure 1. Reflectivity of the corotation barrier in a Keplerian disc. The disc background profiles are given by Q, oc r~ 3 / 2 , Lg = 
Lg c (r/r c ) — 1 , iV 2 = N^(r/r c )~ 3 / 2 , and the sound speed given by c s /(r c O c ) = 0.1. The horizontal axis gives v + — i [see eqs. (132 I t - 
1331 1 . The parameter v is varied by changing the density profile £ oc (r/r c )~ p . The reflectivity for barotropic fluids (Tsang & Lai 2008) 
is recovered for iV 2 = and Lg * = 0, shown as the solid line. 



4.1 Background Disk Structure 

As an illustration of the global p-mode calculation, we consider the standard a-disc model. For the inner region of the disc we 
are most concerned with, radiation pressure dominates gas pressure, and the opacity is primarily due to electron scattering. 
However it is well known that with the standard viscosity prescription for the viscous stress tensor a r( p — — aPtotai, this inner 
disc solution is thermally unstable (Shakura & Sunyaev 1976). We therefore also consider a slightly modified disk model where 
o>0 = — aP S ns is adopted - this disc solution is thermally stable in the inner region (Lightman 1974). Consistent with our 
perturbation analysis, we use the Paczynski-Wiita potential $ = —GM/(r — 2r 9 ) (where r g — GM/c 2 ) to mimic the general 
relativistic effect (Pacyznski & Wiita 1980). 

With the viscosity prescription o r ^ = — aPtotai the relevant background disk profiles are 

E(r) = (6.95gcm- a )MioMi l a- l f^y / Vl-^ 3 ( r i-5trV- 1 (40a) 



P(r) = (3.24 x 10 21 erg cm-^Mfo^ga- 1 (s.) 3/2 (i _ Hla) 1 J (40b) 



H(r) = (1.06 x 10 6 cm)A-/i 8 (l - £2.) ' (l - J . (40c) 

where M 10 = M/(10M Q ), M = (10 18 gs _1 )Mi 8 is the accretion rate, H(r) is the vertical scale height, and J = 1 — l a /l(r) 
with l(r) = fir 2 . The constant l a specifies the specific angular momentum absorbed at the inner edge of the disc per unit 
accreting mass; one typically expects l a ^ l(r[ n ) [l a = l(r[ n ) is the so-called zero-torque condition]. 
The viscosity prescription with a r ^ = — aP gas yields 



E(r) = (7.02xl0 4 gcm- 2 )Mr o 5 M>-5 (£) (l-H^ 5 1 - |f 5 Js , (41a) 



1 1 3 



'., - I _\ ' / 1 _ zLz. \ I I 

3r 

12 l.\ 4 , x -li s -21 



P(r) = (3.27 X 10 25 ergcm- 2 )M 10 5 M 18 5 a"5 (ji) 6 (l - 2i) 5 {l-^JT , (41b) 
£T(r) = (1.06 x 10 6 cm)Mi8 ( 1 - 2?-) ~" f 1 - ^ V . (41c) 



In both cases we obtain the 2-dimensional adiabatic sound speed by c s = TP/E, giving 



^- = 0.72r 1 / 2 M^ 1 M ls (s-) J. (42) 
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Figure 2. The background disc solutions for a-disc models. The depicted profiles are: the vortensity ( = k 2 /(2T,Q) and the modified 
vortensity f e ff = (/S 2 ^ r [both in arbitrary units; see eqs. [1 2 2 1 and i'24l) j: the squared radial Brunt- Vaisala frequency N% , the inverse 
entropy length-scale Lg 1 [see eqs. l|12ll - Ill4jl ] . and the sound speed c s . The adiabatic index is assumed to be T = 1.4. The solid and 
dashed lines denote different angular momentum eigenvalues l . The thick lines show the profiles for the cr r = — aP gaB prescription, 
while the thin lines show the profiles for the o r § = — oPtotal prescription. 



For simplicity, we adopt Y = 1.4 in our calculations (using somewhat different values do not affect our results in Section 4.2). 
Fig- [21 depicts the disc background profiles important for our p-mode calculations. 



4.2 Growing Eigenmodes 

In addition to the outgoing boundary condition (|36[) at some r out 3> ?"ol, it is necessary to impose an appropriate inner 
boundary condition (at n n = nsco) in order to calculate the global p-modes trapped between n n and nL- Unfortunately 
this inner boundary condition is uncertain: the large radial velocity of the transonic flow around nsco leads to energy loss 
of the wave, while the sharp density gradient at nsco provides a partially reflecting inner boundary (see Lai & Tsang 2009); 
in real black-hole accretion flows, a large magnetic flux accumulation inside nsco can make the inner disc edge an even 
better reflector for waves. Here, to focus on the role of the corotational instability, we adopt the free boundary condition (zero 
Lagrangian pressure perturbation) at r- m = nsco, i-e. 

AP= (V + r^MM =0 . (43) 

V dr u Jr ln 

With eq. (|36|) and eq. (|43[) we employ the standard shooting method (Press et. al 1995) with eqs. (|6a[) and (|6b|) . to 
solve for the eigenvalue uj — uj r + iwj. Table 1 summarizes the results for different background disc parameters, and example 
wavefunctions for the zero and one-node eigenmodes are shown in Fig. [3] We see that the trapping region extends from the 
inner edge of the disc at nsco to the inner Lindblad resonance rrL- The wave is evanescent in the corotation barrier region 
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r/r g r/r g 

Figure 3. Eigcnunctions of disc p-modes with m = 3. The disc model parameters are Mig = 3, and l /l(ri n ) = 0.95, with the viscous 
stress tensor given by <j t m = — «P gas . The left panel shows the propagation diagram and eigcnfunctions for the zero- node (in the trapping 
region) mode, with eigenfrequency of o;/f2; n = 2.18 + i3.0 X 10 — 3 . The right plot shows the propagation diagram and eigcnfunctions the 
single-node mode with o;/Q; n = 1.82 + «2.3 X 10 — 4 . For the eigcnfunctions, the solid lines denote the real parts, while the dashed lines 
denote the imaginary parts. 



between riL and rrjL, and tunnels out to the propagation region (r > tol). In the following we will focus on the 0-node modes 
since they have growth rates much larger than the 1-node modes. 

From Table 1 we see that for a given disc model, the (real) mode pattern freqency u} r /m increases only slightly as m 
increases, while the growth rate more rapidly increases with increasing m. In particular, the m = 1 mode has a much smaller 
growth rate than the higher-m modes. These features can be easily understood by examining the propagation diagram (Fig. [4]). 
For small m, the wave trapping region between nsco and til is slimmer, thus to "contain" the same number of wavelengths 
in the trapping region, the pattern frequency must be lower. On the other hand, the wider corotation barrier for small m 
implies that only a small amount of wave energy can tunnel through the barrier, giving rise to smaller corotational wave 
absorption and a slow mode growth rate. As can be seen from Fig. [4l the difference in the width of the evanescent regions 
is greatest between the m — 1 and m = 2; the lower pattern frequency of lower m modes also helps to widen the evanescent 
region. These explain why the m = 1 mode has such a small growth rate compared to the other modes. 

Table 1 shows that for the a-disc models considered, the mode frequency decreases slightly as M increases (while keeping 
the other disc parameters fixed). This results from the increase of the disc sound speed c s [see eq. @2)]. Table 1 also shows that 
the m = 2 and m = 3 modes in each disc model have roughly 2:3 commensurate frequencies, ranging from uJm=3/^m=2 — 1-57 
to 1.69 for the disc background models considered. This has implications for the observations of high-frequency QPOs (see 
section 5). 

Fig. [5] shows the propagation diagram of m = 2 modes and the effective vortensity gradient profiles for both viscosity-law 
(cr r = — Q-Ptotai vs a r< j> = — Q-Pgas) background disc models, with Mig = 3, l /l(ri n ) = 0.95 and T — 1.4. The zero-node modes 
for both models occur at Re(w) ~ 1.34f2j n , which has a positive effective vortensity slope at the corotation radius. However 
only the disc model with the a rt f, — — aP to tai prescription has a growing 1-node mode at Re(u) ~ 1.04f2i n (see Table 1). For 
the disc model with the oy^ = — aP gas prescription, such a (real) frequency would give negative effective vortensity gradient 
at r c , thus the corotation acts to damp the mode (see the inset of Fig. [5]). 

Overall, our numerical calculation of the global disc p-modes is in agreement with our analysis given in Section 3, i.e., 
wave absorption at the corotation resonance gives rise to growing p-modes when the gradient of the effective vortensity is 
positive at corotation. 
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Table 1. Overstable p-modc frequencies for various disc models. 











Mode Eigenfrequencies (u)/Qi n ) for 


O" r — 


— aPtotal 












mode 


P = 


0.9, Mis = 1 


P = 


0.9, Mis = 3 


P = 


0.95, Mis = 


1 


P = 


0.95, Mis = 


3 


m = 


1 


0- node 

1- node 


0.59 - 
0.37 


4- i2.2 X 10" 7 
f i2.4 X 10~ 9 


0.39 
0.17 


+ i4.0 x 10~ 6 
+ i3.1 x 10~ 9 


0.68 


4- i2.2 x 10" 


-8 


0.52 


4- il.4 x 10" 


-7 


m = 


2 


0- node 

1- nodc 


1.46 
1.18 - 


f i4.0 x 10~ 4 
4- i4.2 x 10~ 6 


1.13 
0.80 


+ il.3 x 10~ 3 
+ i5.6 x 10" 5 


1.60 
1.37 


4- il.2 x 10" 
+ i5.5 x 10" 


-4 
-7 


1.35 
1.04 


4- i6.2 x 10" 
+ i2.5 x 10" 


-4 
-6 


m = 


3 


0- node 

1- node 


2.34 
2.01 


4- «2.6 x 10~ 3 
4- tl.a x 10~ 4 


1.91 
1.48 


+ i5.5 x 10~ 3 
+ i7.7 x 10~ 4 


2.52 
2.25 


4- il.2 x 10" 
+ i4.4 x 10" 


-3 
-5 


2.20 
1.83 


+ i3.4 x 10" 
+ i2.6 x 10" 


-3 
-4 


m = 


4 


0- node 

1- nodc 


3.23 
2.86 


4- i6.0 x lO- 3 
4- i9.0 x lO" 4 


2.71 
2.20 


+ il.O x 10~ 2 
+ i2.4 x 10~ 3 


3.44 
3.15 


4- i3.2 x 10" 
4- i3.2 x 10" 


-3 

-4 


3.06 
2.63 


4- i7.3 x 10" 
+ il.2 x 10" 


-3 
-3 










Mode Eigenfrequencies (ui/Q[ n ) for 


= 














mode 


P = 


0.9, Mis = 1 


P = 


0.9, Mis = 3 


P = 


0.95, Mis = 


1 


P = 


0.95, Mis = 


3 


m = 


1 


0- node 

1- nodc 


0.58 
0.37 


4- il.l x 10" 7 
4- i9.6 x lO" 10 


0.37 
0.16 


+ i2.3 x 10~ 6 
+ il.l x 10~ 9 


0.67 


+ i2.3 x 10" 


-8 




- 




m = 


2 


0- nodc 

1- node 


1.44 
1.18 - 


4- »3.6 x 10~ 4 
4- i4.2 x 10~ 6 


1.11 
0.79 


4- il.l x 10~ 3 
+ i6.0 x 10~ 5 


1.59 
1.36 


4- il.O x 10" 
+ i4.3 x 10" 


-4 
-7 


1.33 


+ i4.9 x 10" 


-4 


m = 


3 


0- node 

1- node 


2.32 
2.00 


4- i2.4 x 10~ 3 
4- il.9 x 10~ 4 


1.88 
1.47 


4- i5.1 x 10~ 3 
+ iS.l x 10~ 4 


2.51 
2.25 


4- il.l x 10" 
+ i4.4 x 10" 


-3 
-5 


2.18 
1.82 


4- i3.0 x 10" 
+ i2.3 x 10" 


-3 
-4 


m = 


1 


0- node 

1- nodc 


3.21 
2.85 - 


4- ib.7 x 10" 3 
4- i9.3 x 10~ 4 


2.68 
2.19 


+ il.O x 10~ 2 
+ i2.5 x 10~ 3 


3.44 
3.14 


+ »3.1 x 10" 
4- i3.2 x 10" 


-3 
-4 


3.04 
2.62 


4- i6.8 x 10" 
+ il.2 x 10" 


-3 
-3 



A dash indicates that no growing eigenmode could be found, and that the mode is damped. Here the 
parameter ft = l /l( r in) determines the inner torque condition for the background disc. The mode 
frequencies are independent of the value of the viscosity parameter, a, used. 



5 DISCUSSION 

We have studied the effect of corotation resonance on the adiabatic diskoseismic p-modes (inertial-acoustic oscillations) of 
non-barotropic accretion flows around black holes. Our WKB analysis of the reflectivity of the corotation barrier (Section 3), 
as well as our numerical calculation of the global disc p-modes (Section 4), show that the corotational wave absorption can 
be significantly modified by the non-barotropic effect. In particular, we have showed that super-reflection is achieved when 
[see eq. J33J and Fig. Q] 

v + li - \ > 0, (44) 
where v, /i are defined by eqs. (I20|l - (|21| l. For thin discs, |/x — 1/2| 1 and this condition is simply v > 0, or 

^ lnC -dp7* : = * ln ^J >0 ' (45) 

where C, — k 2 / (2fiE) is the disc vortesnity, A^r is the radial Brunt- Vasala, frequency, and the first equality holds only when the 
adiabatic index F = constant (in which case S = P/£ r ). Thus, in the presence of a reflecting (or partially reflecting) boundary 
at the disc inner edge (r = n n ), the non-axisymmetric p-modes trapped between n n and the inner Lindblad resonance radius 
rm can grow due to corotational wave absorption, when the effective vortensity, £ c fj = (/S 2 ^ r , has a positive slope at the 
corotation radius. As in the case of barotropic discs (Tsang & Lai 2008; Lai & Tsang 2009), the general relativistic effect, 
where k 2 is non-monotonic and becomes smaller as r decreases toward nsco, plays a crucial role in the instability. Now for 
non-barotropic discs, the entropy gradident also plays an important role (cf. Lovelace et al. 1999). 

Our calculations of the global p-modes for various disc models (see Table 1) indicate that the m = 3 and m — 2 modes 
(of lowest radial order) have frequency ratio in the range of 1.57-1.69, similar to the approximate 3:2 ratio as observed in 
high-frequency QPOs of black- hole X-ray binaries (Remillard & McClintock 2006). The growth rates for these modes are 
significantly higher than for the corresponding barotropic case (see Lai & Tsang 2009), due to the effect of the entropy in 
the effective vortensity. Although higher-m modes may grow slightly faster, they would be less likely to be observed due to 
averaging out of the luminosity variation over the observable emitting area. The m = 1 mode is found to have a significantly 
smaller growth rate than the m > 1 modes. 

For our simple Q-disc models, we find that the p-mode frequencies decrease slightly (by about 10-20%) as the mass 
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1.2 




Figure 4. The propagation diagram for disc p-modes with various values of m. The solid curve shows the disc rotation profile Q, 
while the various dashed curves show Q + re/m (above the f2 curve) and Q — n/m (below the fi curve) . The three horizontal lines show 
the representative values of the mode frequency [in units of mQi n , where Q- ln = f2(n n )] for m = 1,2,3 (from bottom to top). The 
corotation resonance is determined by uj/m = f2, the inner Lindblad resonance bj ui/m = Q — n/ra and the outer Lindblad resonance by 
Ld/m = Q + n/m. A mode is trapped between r- ln = nsco an d nL, an< ^ ' s evanescent between hl and roL (the horizontal dotted lines). 



accretion rate increases by a factor of 3. Observationally, it is known that high-frequency QPOs are observed only when 
the X-ray binary systems reside in the so-called steep power-law spectral state (also called "very high state"), which may 
corresponds to a very specific range of accretion rates. It is unclear whether our result is consistent with the observed trend in 
high-frequency QPOs (e.g., Remillard et al. 2002; Remillard & McClintock 2006). Clearly, more sensitive observations (e.g., 
with future X-ray timing missions; see Barret et al. 2008, Tomsick et al. 2009) would be useful to determine this trend and 
to search for the possible m = 4 mode and the frequency ratios. 

Finally, it should be noted that our calculations of global disc modes are still based on rather crude models. The a- 
discs are phenomenological models, and our results (especially the mode growth rates) depend sensitively on the inner disc 
boundary conditions (both the lo parameter for the background disc and the reflecting boundary condition for the waves). 
Other potentially important effects (such as turbulence) have not been taken into account. Thus we should treat our specific 
results (such as those presented in Table 1) only as a demonstration of the basic physical principles, and any comparison with 
the observations at this point should be taken in this spirit. 
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Figure 5. The propagation diagram for m = 2 disc p-modes and the derivative of the effective disc vortensity profiles. In the lower panel, 
the £ e ff profiles are shown for disc models with the viscosity prescription a Tl k = — c*P gas (solid curve) and <r r = — aPtotal (dashed curve), 
and the other disc parameters are Mis = 3, i D /'( r in) = 0.95. In the upper panel, the two horizontal lines give the real eigenfrequencies 
of the 0-node mode, Re(w) ~ 1.34f2; n (for both disc models) and the 1-node mode, Re(aj) ~ 1.04f2; n (for the oy^ = -aP tota i disc model 
only). Note that growing modes can be found only if dln^ a g /dr > at the corotation, thus the Re(a;) ~ 1.44f2i n mode exists for both 
disc models, while the Rc(w) ~ 1.04f2; n mode exists only in one of the disc models. The inset of the lower panel shows a magnified view 
of the derivative of the effective vortensities at the corotation point of the Re(w) ~ 1.04Q; n mode. 
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